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HESSIANS OF SPECTRAL ZETA FUNCTIONS 


K. Okikiolu 


Abstract. Let M be a compact manifold without boundary. Associated to a metric g on M there are various 
Laplace operators, for example the De Rham Laplacian onp-forms and the conformal Laplacian on functions. 
For a general geometric differential operator of Laplace type with eigenvalues 0 < Ai < A 2 < • . . we consider 
the spectral zeta function Z{s) = 7^0 The modified zeta function Z{s) = r(s)Z(s)/r(s — n/2) is 

an entire function of s. For a fixed value of s we calculate the hessian of Z{s) with respect to the metric and 
show that it is given by a pseudodifferential operator Ts = Us + Vs where Us is polyhomogeneous of degree 
n — 2s and Vs is polyhomogeneous of degree 2. The operators Us/T'{n/2 + 1 — s) and T4/r(n/2 + 1 — s) 
are entire in s. The symbol expansion of Us is computable from the symbol of the Laplacian. Our analysis 
extends to describe the hessian of {d/ds)^Z[s) for any value of k. 


1. Introduction. 

Let M be an n-dimensional smooth compact manifold without boundary. The tensor powers of M"' 
are TP’'^{W^) = (g) (gj'^'R’^*. Suppose U is a vector subspace of TP’'^(R"') which is invariant under 

the action of the general linear group GL(R^). Writing V for the principal frame bundle of M, we 
obtain a general tensor bundle E over M by forming the associated bundle 

(1.1) E = V 

For example we obtain in this way T^’'^(M), the bundle of type-(p, ( 7 ) tensors on M, and S'^M, the 
bundle of symmetric (0, 2)-tensors. Let End(E) denote the bundle whose fibers are the endomorphisms 
of the fibers of E, so End(E) can be identified with E ® E* . Let N denote the rank of E, and let 
C°°{E) denote the smooth sections of E. 

For a Riemannian metric g on M, let dV denote the canonical volume element in the metric g, and 
V the total volume of M. The metric g naturally gives a metric (, )g on E which gives an inner product 
((,))g on C-(E), 

{{9,co))g = [ {e,u;)gdV, e,u;eC°°{E). 

J M 

The metric g also gives rise to the Levi-Civita connection on M with curvature tensor R, and the 
connection on E 

V = Vg:C°^{T'^M®E) 

which is invariant under scaling, that is for c > 0 and uj G T^{M (g) E) we have VcgUJ = VgLO. 
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Now suppose that for each Riemannian metric g we have a differential operator F = Fg : C°°{E) 
C^{E) which is self-adjoint with respect to ((, ))g, non-negative, scale invariant, and geometric of 
Laplace type, meaning that 


(1.2) Fco = [^V\ei,ei,La) + ^ R{ei,ej,ek,ei)A{e\9^,9^9\u;)] + Beo 

i,j,k,£=l 


2=1 


for CO G C°°{E), where R is a constant, ci,..., e^, is a pointwise gf-orthonormal base for TM with dual 
coframe 9^, 9^, and A G C°°{T^’^{M) 0 End(£')) is obtained from some element 


which is invariant under the natural action of Indeed, if O is the gf-orthonormal frame bundle 

of M, then 

T^'^{M)®E®E* = O Xo(K-) ®U*) 

and A naturally gives rise to a section A = Ag oi 0 E 0 E*. Examples of F include the de 

Rham Laplacian on forms and the conformal Laplacian. 

Write for the real part of s G C. Let Ai < A 2 < • • • be the eigenvalues of F and set 

CX5 

Z(,) = ^A-‘, 
i=i 


Then Z{s) converges for iRs > 
simple poles at 

s 


nl2 and extends analytically to a meromorphic function with possible 

f n/2, n/2 — 1 ,... 1, for n even, 

( n/2, n/2 — 1,. for n odd. 


In fact, if F is invertible, then 


Zis) : = 


T{s)Z{s) 


r(s — n/2) 

is entire. In general, when F has a kernel of stable dimension d, then 


(1.3) 


Zg{s) + 


d 

r(s — n/2) s 


is entire. 


and so for two metrics g and g, Zg{s) — Zg{s) is entire. 

Remark. (1.3) does not hold for general poly homogeneous pseudodifferential operators of stable dimen¬ 
sional kernel in odd dimensions, since Z(s) may not vanish at the negative integers. See [Okl] 1.2 for 
further details. 

With s fixed, differentiating Z(s) twice with respect to the metric one obtains the Hessian of Z(s) 
which is a symmetric bilinear form on the space C°°{S‘^M). At the metric g, 


Hess2^(s)(/i, h) 


da'^ 


Zg.\-ah{s). 

(T = 0 


Set N+ = {1,2 ,...}. 
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Theorem 1. For s G C, there exists a unique symmetric pseudodifferential operator Tg = Ts{F) : 

such that 

(1.4) RessZ{s){h,h) = {{h,Tgh))g. 

The operator Tg is analytic in s. For s ^ nl2 + N"*", there exist poly homogeneous pseudodifferential 
operators Ug and Vg of degrees n — 2s and 2 respectively such that Tg = Ug + Vs. The operators Ug and 
Vg are meromorphic in s with simple poles located in n/2-\- N+. (The poles of Ug and Vg cancel in the 
sum Tg = Ug + Vg, but the symbol expansion of Tg for s G n/2 + N+ involves logarithmic terms.) For 
general s, the poly homogeneous symbol expansion of Ug is computable from the complete symbol of the 
operator F. In particular, there is a simple algorithm to compute the term Ug of homogeneity n — 2s. 
Furthermore, we can differentiate in s to obtain 

(1.5) Ress{d/ds)'^Z{s)ih,h) = {{h,{d/ds)^Tgh))g, 

and the principal symbol of {d/ds)^Ug is equal to the leading order term of {d/ds)^Ug (provided this does 
not vanish identically). 

Theorem 2. The symbol Ug from Theorem 1, can be computed as follows. Write F' = {d/da)\a=QF{g+ 
ah) and let x G M. Take coordinates on M which are orthonormal at the point x, and take a local 
trivialization of E on a neighborhood of x to obtain coordinates for E. Suppose that in these coordinates, 
the operator F' is given at the point x by 

( 1 - 6 ) ^ ai, 

where a and (3 are multiindices and A)() ^ is an N x N real valued matrix. Set 

S = s — n/2, 


write I for the identity operator on E,^, and set 


(1.7) 


C(s) 


. 1 y/2 

r (-25 + 2 )‘ 


Then at (x,^) G T*M, the value of Ug{x,f/) G End(5^M)a; is given by 


(1.8) {ug{x,^)h)ij = ^ Ug{d°‘,d^,d'^,d\x,^)tiace{A)/^^{x)A'(^^^{x))hki, 

H +\ 0\=2 
\j\ + \S\=2 

k,e 

where the terms Us{d°^,d^,d^,d^,x,f/) are given as follows: 

(a) Ug{djdk, I, dpdq, I, x,f,) = 4(45^ - . 


(b) 


ngi^dj, dk, dpdq, I, x,(/) 


- 2 ( 45 ^ - 
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(c) Us{dj, dk, dp, d„ x,0 = (452 + 25 - - (25 - 

(d) Us{I,d,dk,dpd„I,x,0 = (452-25K,a^p^,|er-2^-" + (25-l)5,,U,|^r-2^-2. 

(e) Us{l,d,dk,dp,d„x,0 = -{‘^S^ + s 

+ (5 - l/2){-5,kU, + 

(f) Us{I, djdk, I, dpdq, x,C) = 

{s^+s)i,^k^pUT~^^~^ + \{s-i)5,kU,W~''^~^ + \{s-i)^,ikSpq\^r‘^^-^ 

~ l^i^jp^k^q ^kq^j^p + ^jq^k^p + ^kp^j^q)\^\^ 2s 2 _|_ -^{^SjkSpq + 5jp5kq + 5jq5kp)\^\^ • 

We refer to the operator Tg in Theorem 1, as the operator corresponding to Hess^(s). At a metric 
g which is critical for Tg, the eigenvalues of Tg can determine whether gf is a local maximum, minimum 
or saddle point. Computation of the principal symbol of Tg gives at least some information about the 
type of critical point, for example if it changes sign, one can say for sure that is a saddle point, while 
if it is positive, one can generally determine that g has finite index. The operators Ug, 14 are only well 
defined up to the addition (respectively subtraction) of a smoothing operator which is analytic in s. 
They are however canonically defined on the level of symbol expansions because in general they have 
degrees which differ by a non-integer. It should be noted that while the symbol expansion of Ug is 
locally computable, the expansion of 14 generally involves global quantities. Theorem 1 sheds light on 
the question of [GS] concerning how many derivatives of an invariant must be computed to obtain a 
local expression. We see that for low values of s, it takes two derivatives of a zeta invariant to obtain 
an expression whose “leading order” is local. Indeed, in general Ug is the principal symbol of Tg when 
iRs < n/2 — 1 and this is local, while Vg is the principal symbol when iHs > n/2 — 1. Of particular 
interest in this discussion is the zeta-regularized determinant of F, which is defined to be e~^ The 
quantities u(0) and {d/ds)\g-oUs usually control the type of critical points possible for the determinant. 
The picture may become complicated at values of s for which Ug or Vg vanishes. At these values, the 
principal symbol of Tg is given by the highest degree non-vanishing term in the symbol of Ug or 14 , or if 
these terms have the same degree it is given by the sum. For example, in [Ok2] we will actually compute 
the principal symbol of 14 for the family of operators in (1.9). We will see cases when 14 has order zero, 
and Ug is the principal symbol of Tg for JRs < n/2. We also remark that one can formulate a slightly 
more abstract but more basic version of Theorem 1, which does not mention metrics or geometric 
operators. Indeed, fix a tensor bundle over M, and consider the space of elliptic, non-negative second 
order differential operators with isotropic principal symbols. Now consider the space of coefficients of 
such operators and compute the Hessian of the zeta function with respect to variation of the coefficients. 
Then this bilinear form on the space of coefficients is given by a pseudodifferential operator, in a similar 
fashion to Theorem 1. Theorem 1 then follows from the more abstract formulation by expressing the 
coefficients of a Laplace type operator in terms of the metric. 

In this paper we prove Theorems 1 and 2. Although similar results hold for geometric elliptic 
operators of more general form, we will not deal with greater generality here. For us the main goal is 
to understand the behavior of zeta functions for specific geometric operators. In two further papers 
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[Ok2] and [OW], we carry out specific calculations. In [Ok2], we explicitly compute Us for the family 
of Laplacians 

(1.9) F = + ci^5) + C 2 I, 

where A is the (positive) Laplace-Beltrami operator, S is the scalar curvature, I is the identity operator, 
IJ, = (n — 2)/(4(n — 1)), and ci,C 2 are arbitrary constants which render F non-negative. This family 
contains both the Laplace-Beltrami operator and the conformal Laplacian (ci = 1, C 2 = 0), and we are 
in fact able to compute the principal symbol of I 4 which is non-local. We go on to make conclusions 
about the critical points of Z{s) for all values of s. For large values of s, Z(s) should be dominated by 
the lowest eigenvalues, and one sees consequenses of this in the results. In joint work with C. Wang, 
[OW], we compute Ug for the de Rham Laplacian. Theorem 2 (a)-(f) gives the symbol Ug needed for the 
applications [Ok2] and [OW]. The general formula for Ug is given in Lemma 5.1. We do not give here 
explicit formulas for the complete symbol expansions of Ug and 14 , although the symbol expansion of 
Ug can be obtained from (5.10) or Lemmas 6.1, 6.5, and the symbol expansion of I 4 has contributions 
from three terms, firstly (4.5), secondly (4.16) whose symbol expansion can be computed from (4.18), 
and thirdly (4.8) whose symbol expansion is computed in Lemmas 6.1, 6.5. 

This work is motivated by a growing list of results [OPSl], [BCY], [CY], [Br], [Ri], [Okl], [Mol], 
[Mo2], where the authors take F to be the conformal Laplacian or Laplace-Beltrami operator, and 
identify local or global maximal or minimal metrics for the functional logdetT or Z(s), either within a 
conformal class, or among all smooth metrics. Papers [OPS2], [OPS3] and [CGY] contain applications 
of some of those results. It emerges that in some situations only maxima occur while in others only 
minima occur. In [Ok2] and [OW] we will address the questions of when it is possible for maxima to 
occur, when it is possible for minima to occur, when neither can happen, and when one might hope to 
find a convexity argument to prove global results. 

Acknowledgements. The author would like to thank the referees for their careful and extensive 
comments on earlier versions of this manuscript. In addition, the author would also like to thank the 
following institutions for support: ESI, IAS, IHES, MSRI, Princeton University. 

We will now indicate briefly what is involved in the proof of Theorems 1 and 2. One can use the 
parametrix for the heat kernel or the resolvant of the Laplacian to prove the result. Technically the 
work involved is comparable for the two methods, and each method has its advantages. The heat kernel 
method has the advantage of yeilding completely explicit formulas for the symbol expansion of Ug, and 
so we take this approach. We fix a metric g, and take a smooth one-parameter deformation g{<T) with 
ff(0) = g- Using the Melin transform to write in terms of the heat kernel, we can differentiate 

to compute the second derivative of Zg{s) with respect to a at a = 0, see Lemma 3.1. The main term 
we need to study is the one giving rise to the operator Ug, which is the last term on the right hand side 
of (3.1). Writing F' = {dFg(^^'^/da){U), it equals 

(1.10) ——^^ trace f f (u + v)'‘‘F' F' dudv. 

T{s-n/2) JJu+v<i 

Before proceeding, we give a list of the steps in the proof to help the reader navigate. The steps in this 
guide will usually be referenced when they occur in the proofs that follow. 

Table 1.1: Steps of the Proof. This table lists the steps which are needed to go from the expression 
(1-10) to the symbol expansions of the operators Ug and 14. They are not listed here in the order they 
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occur in the proof, and some of the steps need to he repeated more than once in the proof Steps 1-8 
are entirely trivial while Steps 9-16 represent the analytic details, with Steps 13-16 being the more 
significant or subtle. On the other hand, one has to pay a little attention to steps 1, 8 and 9, because 
they change the nature of the expression being considered. 

Step 1 ; The steps just earried out produeed a sum of terms. Consider each term separately. 
Step 2 ; At a previous stage, the expression was decomposed as a sum and each term was 
considered separately. Reverse this by summing up the eontributions from each term. 

Step 3 : Introduee a partition of unity to reduce to coordinate eharts over which E is trivial. 
Step 4 : Work in geodesic normal coordinates. 

Step 5 : Write h = ff^(O) in eomponents. 

Step 6: Write F' in components. 

Step 7: The expression has the form “ff ((u -h v)^ - dudv”. Restriet attention to the term 

“ -” in the integrand. 

Step 8; Reverse Step 8 by reinstating the faetor {u + v)^ and the uv integration. 

Step 9; Identify an operator T (by its kernel) so that the expression has the form {{h,Th)). 
Step 10; Write the traee of an operator as the integral of the kernel over the diagonal. 

Step 11; Make the standard small time expansion of the heat kernel. 

Step 12; Make a Taylor expansion of eaeh term in a given variable. 

Step 13; Colleet together terms of the same “order”. 

Step 14; Exchange uv integration in the expression with another limiting proeess sueh as taking 
the operator trace, a Fourier transform or integration in other variables. 

Step 15; Compute the symbol of the operator from its kernel by taking the Fourier transform. 

Step 16; The symbol has the form ff((u-l-v)^ - dudv”. Carry out the uv integration to 

see the asymptotie behavior in the eotangent variable 

The simplified scalar case; We wish to write (1-10) in terms of /i = {dg/da){0). Let us first pretend 
for the sake of clarity that F' is just multiplication by the smooth function cj) = hijA^^, where A^^ 
is a smooth tensor field. In doing so we will be ignoring some technicalities which arise when yl is a 
differential operator, but let us concentrate here on the main point. This eliminates the need for Steps 
5, 6. We now wish to compute 


( 1 . 11 ) 


r(s — n/2) 


trace 


{u-IvYfie e dudv. 


u-\-v<l 


Steps 14, 10: Writing K{e *^,x,y) for the integral kernel of e this is equal to 
1 


( 1 . 12 ) 


{u + vf 


4>{x) K{e ,x,y) (j){y) K{e ,y,x) dV{x)dV{y) dudv. 


—vF 


r(s - n/2) JJu+v<i J JMxM 

Steps 14, 9: We see that this equals 
(1.13) 

where Tg is the operator on functions with 


^,TgC 


(1.14) 


ifs(j)){x) = [ Kg{x,y)(j){y)dV{y), 
J M 
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where 


(1.15) Ks{x,y) = — — ff (u + vY K{e '^^,x,y)K{e ,y,x) dudv. 

r(s - n/2) JJu+v<i 

To express this operator as a psuedodifferential operator, we compute its symbol by taking the Fourier 
transform. Here, we focus just on the leading order term in this symbol. 

Steps 3, 1, 4: Fix a point x G M and work in normal coordinates w of y about x. 

Steps 7, 11, 13: For this discussion, let us assume V = 1, and restrict attention to the leading order 
term of the heat kernel K{e~^^,x, y), which a constant multiple of, Approximately then 

we have 

(1.16) K{e-^^,x,y) K{e-^^,y,x) ^ . 

Step 15: Since '^J~g(x) = 1, up to leading order, the symbol of the right hand side of (1.16) is obtained 
by taking the Fourier transform in w. This yields a constant multiple of, 

(1.17) (^i + ^)W2g-|«PW(«+«). 

Steps 14, 8 : From this, we see that at x, the symbol of T^, whose kernel is given by (1.15), is 
approximately 


Step 16: Substituting t = u/{u + v) and T = |^P(ii + v), this is 


(1.19) 


F(s — n/2) JT<\i\‘^ Jt=o 


j.s-n/2+lg-Tr(l-r) 


If we replace the interval of integration by 0 < T < oo, we get 


i-2s-4 /-CO ^1 


( 1 . 20 ) 


F(s — n/2) 


pOO pi 

Jt=0 Jt=C 


j.s-n/2 + l^-Tr{l-T) ^ 


F(n/2 - s - 1) F(n/2 - s + 1) 1^/ 
F(n — 2s — 2) 


This gives the principal symbol of Us- The reason it has degree n — 2s — 4 instead of n — 2s as stated 
in Theorem 1, is because we have assumed that F' is just multiplication by the function cj) = hijA^Y 
whereas in general it may involve derivatives up to order 2. Since F' occurs twice in (1.10), this accounts 
for the four missing derivatives. To compute (1.19) we need to subtract from (1.20) the part of the 
integral on the left hand side over the interval < T < oo. But by steepest descent, there is an 
expansion as T ^ oo of the form 


( 1 . 21 ) 


-Tr(l-r)^^ ~ 2r-^ + CoT-^ + C^T-^ + 


Hence as |.^| ^ oo. 


( 1 . 22 ) 


T{s — n/2) Jt>\p\^ 


rjns — n/2-\-l 


-Tr(l-r) 


-1 ( 2|C|-^ C2|^|-^ 

F(s — n/2) \s — n/2 + 1 s —n/2 s — n/2 — 1 
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This is a symbol of order —2. The poles, which occur at s = n/2 + 1, n/2 + 2,..., cancel with the poles 
of (1.20) and contribute logarithmic terms to the symbol of Tg at these values of s. 

The general case: Now we modify this basic calculation to obtain Theorem 2 on the tensor bundle 
E. 

Step 5: Suppose that h is supported on a small neighborhood of a point x on which we take 
coordinates. 

Step 6: Write F' in components on Q as 


so that (1.10) becomes 
1 


(1-23) 




r(s — n/ 2 ) 




trace / / {u + vYA^l^ {d‘^hij)d^ e {d"^hki)d^ e dudv. 

J J u-\-v<l 


In fact, because we are only interested in the leading order term in the symbol, we can restrict the sum 
in (1.23) to those multiindices with \a + f3\ = \l + d\ = 2. 

Steps 1, 14, 10: We will analyze each term in the sum (1.23) separately. We write 


(1.24) 


1 


trace 


r(s — n/2) 

r(s-n/2) 


u-\-v<l 

u + v) 


{u + v)^AY 0 (d^hiYd^ A'/i {d^hu)d^ dudv 

{d°‘hij){x) E{u,v,x,y) {d'^hki){y) dV{x)dV{y) dudv, 


MxM 


where 

(1.25) E{u,v,x,y) = trace A"lp{x) ,x,y)) A/^l{y) {d^yK{e-'’^,y,x)). 

Steps 14, 9: This has the form 


(1-26) {{h^.Xhhkd), 

where T^lj is the operator on scalar valued functions whose kernel to leading order 


IS 


(1.27) 


r(s - n/2) JJu+v<i 


{u + v)^d)/d/E{u,v,x,y) dudv. 


{We integrated by parts to obtain this. Derivatives falling on the volume factors \/g{x) and \/g{y) 
yield lower order terms.) 

Steps 11, 13: To compute the principal symbol of the kernel in (1.27), we restrict attention the leading 
order term of the heat kernel K{e~^^ ,x,y), which is a constant multiple of, /44 where 

d{x,y) is the distance from x to y. This enables us to replace (1.27) by 


(1.28) 


(-1)1/3+71 trace yl)/_g(x) A^^x) 


r(s — n/2) 


u-\-v<. 


(n + u)* (nu )-^/2 


X 


dudv. 



Here, we are using the fact that to leading order, 

Qa^-d?{x,y)lAu\ ^ f — iV"! 

X \ x=0 ^ ' y ^ 

Step 4 : At this point we fix the point x and work in geodesic normal coordinates on H centered at x. 
(The coordinates used in (1.23) to reduce to a scalar problem may be arbitrary.) 

Steps 14, 15, 16: Taking the Fourier transform of (1.28) in the normal coordinates of y as in (1.17)- 
(1.20), yields the degree n — 2s term in the symbol of T^lj of (1.26). Carrying out the calculation for 
explicit choices of {a,(3,'y,6) gives rise to the quantities in Theorem 2, (a)-(f). 

The rest of the paper is organized as follows. In Section 3, we derive variation formulas for Z(s) to 
write the Hessian of 2{s) in terms of derivatives of F. In Section 4, we use basic estimates on the heat 
kernel to reduce the problem of studying the Hessian to that of studying (1.10) (equivalently (4.19)). In 
order to analyze the operator associated to (1.10), we will need to make a double decomposition of the 
integral kernel, one stage to extract the polyhomogeneous expansion of the symbol of Ug , and the other 
to extract the polyhomogeneous expansion of I 4 . In Section 5, we extend the above sketch to obtain 
the expansion of the symbol of Ug and use this to deduce Theorem 2. However, we leave the rigorous 
estimates on errors to Section 6 where we prove Theorem 1 by carrying out the double decomposition 
in detail. Formula (1.5) then requires no proof because it follows directly from (1.4), the fact that Tg is 
an analytic family of operators, and the standard fact that if the kernel of F has stable dimension then 
Zg{s) is smooth on A4 x C. 

2. Background and conventions. Before commencing with the proof of Theorem 1, we discuss in 
more detail the Hessian of the function Z{s) and we fix the conventions which will be assumed concerning 
polyhomogeneous pseudodifferential operators (TDOs). Let A4 denote the space of Riemannian metrics 
on M. Define the Sobolev space F[^{S‘^M) to be the completion of {S‘^M) in the norm || ||^ defined 
by 

\\h\\r = Y.{{v^h,v^h))y^ 

s<r 

where is a fixed but arbitrary smooth metric. The Sobolev imbedding theorem shows that for r > n/2, 
is contained in the space C{S‘^M) of continuous sections of S'^M. The smooth topology on 
M is the Frechet space topology coming from the collection of all the norms || H^.. For r > n/2, let A4'’ 
denote the closure of the space of smooth metrics on M in the space For a compact subset 

U C C, let Hol(C/) be the space of holomorphic functions on U with the supremum norm. 

Proposition. Suppose m, r G N with r > m + n/2 + 1, and 1/ C {s G C : > m + n/2 + 1 — r}. Then 

the map g Zp^ is an m-times eontinuously differentiable function from A4’’ to Hol(t/). 

The proof is given in [Ok3]. Using it, we obtain the Hessian of {d/ds)^Z{s) as a bilinear form on 
F[^{S‘^M), provided r > n/2 + 3 — min{(0,3f?s)}. This restricts to a bilinear form on C°°{S‘^M). We 
should remark that the Hessian defined above Theorem 1 is computed using the local linear space 
structure of A4, that is using the flat connection on A4. It is possible to take other connections on A4. 
For a sufficiently regular function F : M ^ C and a smooth curve g{cr) in A4, the derivative DF is the 
linear form on C°°{S‘^M) defined by 


DF 



dF{9{(^)) 

da 
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The Hessian of F relative to a connection V on is given at G by 


( 2 . 1 ) 


Hess F 


\da' da) 


d^F{g{a)) 

da'^ 


DF V 


\da' da)) 


We will prove Theorem 1 when the Hessian is defined relative to the flat or intrinsic connection 
on Ai. The intrinsic connection on A1 is the Levi-Civita connection coming from the intrinsic metric 
defined as follows. For u,v e C°^{S‘^M), at the point g G M, 


Then 

( 2 . 2 ) 


^intrinsic 


Vl^,^ 

' da'da 


d^g. 






da"^ 


~ke^9ik^^ , 1 ~k(. dm dgjj 

^ da da 2^ da da 


-kl-pqpdip dgkq 




Now we give a summary of some facts we will need about polyhomogeneous TDOs. For further 
details the reader is referred to [Tr] and [Se]. If [/ C M"" is open and m G M, a function q G C°°{U x 
M"', r^’^(R^)) is an iV X N-matrix valued symbol of order m (and type (0,1 )), if for all multiindices a, (5 
and compact sets K C U, 

md)q{x,0\ < C(l + |^|)"^-|^l, 

where C = C{a, j3, K) and | • | is a norm on r^’^(M'^). The symbol q has order —oo if it has order m 
for all m G M. The symbol q is polyhomogeneous symbol of degree // G C if there exists a sequence of 
functions qj G C'°“(n x G ^ / 0}, j = 0,1, 2,... , which are homogeneous; 

qjix,rf) = rf^-^qj{x,f,), r > 0, 

such that 

q{x,f) ~ qo{x,f) + qi{x,^) + ..., as |C| ^ oo. 

By this we mean that if we introduce a cut-off function f G (^“([O, oo)) supported on [1/2, oo) with 
(/ = 1 on [l,oo), then 


.7-1 

q{x,f) - (pilCl) qj{x,f,) is a symbol of order Ifig — J. 
j=o 


(Where is the real part of ^.) 

Now suppose TT : FI ^ n is a rank N vector bundle and Q : C°°{E) C°°{E) is a linear operator. 

Then Q is a TDO of order m if whenever we can trivialize E over a coordinate neighborhood thus 
identifying it with (U C M”) x then there is a symbol q of order m such that 

(2.3) Qh{x) = —^ [ [ e^'^'^~y'>'^q{x,f)h{y)dydC. 

for X in 17 and h supported on a compact subset of U. The principal symbol is then a section of the 
pull back of the bundle E 0 E* to the cotangent bundle of M. Q is polyhomogeneous of degree fi if 
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and only if for all trivializations of E, the symbol g of Q is polyhomogeneons of degree /x. We remark 
that from this it follows that there is a smooth kernel k on M x M \ {{x, x)} such that 

Qh{x) = / k{x,y)h{y)dy 
Jm 

whenever x is not in the support of h. Conversely, if the operator Q : C°°{E) C^{E) is given by a 

smooth kernel away from the diagonal then it is a 'kDO of order m provided (2.3) holds on coordinate 
charts of small diameter. If the operator Q is given by a kernel which is smooth everywhere, it is called 
a smoothing operator. 

When the manifold M is equipped with a Riemannian metric g, it is convenient to use the density 
dV to define the kernel K{Q,x,y) of the operator Q, so 


Qh{x) = [ K{Q,x,y)h{y)dV{y). 

Jm 

With this convention, the symbol q{x,^) of Q in a given coordinate system is obtained by taking the 
Fourier transform of K{Q,x,y)dV{y) : 


g(x,0 = j K{Q,x,y)x{y)e^^y ^^'^dV{y). 

Here, x is a cut off function equal to one on a neighborhood of x and supported in the coordinate chart. 
The symbol q is well defined up to a symbol of order —oo arising from the choice of x- In particular, 
if Q is polyhomogeneons symbol expansion is well defined at x in the given coordinate system. It will 
generally be convenient to compute the metric-modified symbol 

(2.4) q{x,C) = J K{Q,x,y)x{y)e^^y~^^'^ dy 

rather than the true symbol. The two are of course very simply related. 

Suppose now that we choose a smooth function re : H x H ^ such that for x fixed, w{x, y) defines 
geodesic normal coordinates of y centered at x. Then (x, y) (x, w) defines a smooth diffeomorphism 
from a neighborhood of the diagonal in H x H into H x R”^. It follows from [Ho] that (2.3) can be 
replaced by 

/ / e'^'^j{x,i)h{y{x,w))dwdi, 

whenever there is a set of small diameter in M which contains both the support of h and the point x. 
Here, g is a symbol of order m if we express x in coordinates. 

If Qs : C°°{E) C°°{E) is a family of ^^DOs depending on a parameter s G O C C, then we say 

that Qa is a analytic family of polyhomogeneons ^DOs if Qg has degree /x(s) where y is analytic on O, 
and on any local trivialization of E, Qg has symbol q{s,x,^) which is smooth in (s,x,f) analytic in s, 
and satisfies the estimates 

\dfd^^q{s,x,0\ < C(l + |e|)''^(^)-|^l, 

where ifiy{s) is the real part of yfs) and C, which depends on a, f3, is uniform for x and s in compact 
sets. This implies that the homogeneous terms qj{s,x,f) in the symbol expansion of Qg are analytic in 
s, that the kernel defining Qg away from the diagonal is analytic in s, and that if A : C°°{E) C°°{E) 
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is a positive elliptic second order differential operator, then QgA is a analytic family of bounded 

operators on L‘^{E). 


3. Variation formulas. 

Let g he a metric on M, let h G C°°{S‘^M), and set 

5 ( 0 -) = g + ah. 

We use primes to denote differentiation with respect to a, so for example 


F' = Sf, 


da 




F" = — F; 


da'^ 


■ g{^)- 


In this section we will prove the following Lemma expressing the second variation of Z{s) in terms of 
F, F' and F". 

Lemma 3.1. 

dPZ{s 




trace ( F'UF'F-^ 


+ 


1 


r(s — n/2) 


trace 


(3.1) 


u-hv>l 

+ 


J F (e-‘^ - n) dt 

{u + vfF' - n) F' (e-^^ - U) dudv 

1 


r(s — n/2) 


trace 


{u + v)^F' F' dudv 




We start by expressing the first variation in terms of F and F'. In the case when the kernel of F is 
trivial, 

(3.2) Z{s) = trace F-, 

r(s — n/2) 

and 

DZ{s){h) = -^Zis) = trace F'F-^-^. 

da r(s — n/2) 

When F has a non-trivial kernel whose dimension does not vary with a, let 11 denote the projection 
onto the kernel, and write 

F-^ = (F + n)“* - n, 

so (3.2) still holds. Using 11^ = 11 we get 


and using FII = 0 we get 


n'n + nn' = n', 

F'n = -Fn' 
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so 

nF'n = 0 


and 


trace F'n = 0. 


Hence 


tracen'(F + n) ^ ^ = 2tracen'n(F + H) ^ ^ = 2tracen'n = traceH' = 0, 

and 

DZis)ih) = -£ii±^trace(F + n)'(F + n)-^-i 
i (s — n/2) 


It will be convenient to write this in the form 

(3.3) DZ{s){h) = trace F'L(s) 

where 

Now we compute the second variation of Z(s) in terms of F, F', and F". In the case when the kernel 
of F is trivial, 


-r(^ + i) 

r(s — n/2) 


trace F'{F + H) 


— s — 1 


"^/nl traceF'F * ^ 
r(s — n/2) 


d^Z{s) -r(s + i) 
da^ ~ r(s - n/2) 

Using the Melin transform, we write 


trace F F 


n 7-1—s — 1 


r(s + i) 

^ ’ ■trace F' 


r(s — n/2) 


da 


n —S —1 


1 


r(s +1) Jo 


Fe dt. 


By Duhamel’s principle. 


de 


-tF 


and so 


da 


I 


da 


t 


f ^ e-d-FF du, 

I u=o da 


cr=o do da 

1 


«00 /*CXD 


a=0 ^('5 + 1 ) Jo 

s ^—uF Tpf ^—vF 


r^OO /‘t 


F 6-“-^ F' dudt 


r(s + i) Jo 


{u + vY F' dudv. 


'0 
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Hence 

d^Z{s) 


da^ 


-r(g + i) 

r(s — n/2) 


-I pCXD pOO 

trace + —--—trace / / (u + vY F' e~'^^ F' e~'’^ dudv. 

r(s - n/2) Jv=Q 


In the case when F has non trivial smoothly varying kernel, we have HF = 0 = FH, so we get the 
relations 

FH' = -F'n, n'F = -HFh 


They imply in particular that 

FH'F = 0, 

Since H' = H^H + HH', we also have 


so 


Hence 


We dehne 


HF'H = 0. 

HH'H = 0, 

(F + n)n'(F + n) = fh'h + hh'f = -f'h - hf'. 

H' = -(F + n)-iF'n - nF'(F + n)-^ = -f-^f'h - hf'f-^ 


= { 


f (1-^—i)/(^-l) ^/0,1 

S + 1 M = 1 

0 /U = 0 


Then when /u > 0, 

J 

Starting from the formula 


r »00 /*C >0 


0 ^0 


{u + v)^ e-^ dudv = r(s + l)T,(/r). 
dZ{s) -r(s + i) 


da r(s — n/2) 

and following the previous argument gives 
d^Z{s) _ -r(s + i) 


trace F'(F + H) 


— S—1 


(3.5) 


trace F"(F + H) 


— S —1 


da^ r(s — n/2) 

-I pOO pOC 

+ — - —trace/ / (u + u)"F'(F + H)'dudu. 

r(s-n/2) Jo Jo 


Furthermore, 

pOC pOO 

trace / / {u + v)^F' H' dudv 

Jo Jo 

POO POO 

= -2trace / / {u + vf F' H F' F-^ dudv 

Jo Jo 


>0 Jo 

= -2trace ( F' H F' F“^ 


^OO POO 


'0 -'O 


= -2r(s + 1) trace H F' F"^ 'i’,{F)F', 
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and since 


^-u{F+ii) ^ + e““n, 

we get 

F' F' = F' (e““^ - H) F' - H) 

+ e““ F'n F' + e“"'F'e““^F'n, 

and 

nOO pOO 

trace / / (« + z;)"F' F' dudr; 

^oo ^oo 

= trace / / (w + (6”“^ - H) F' - H) dudr; 

JO Jo 

+ 2r(s + 1) trace IIF''I's(F) F'. 

Moreover 

trace F"n = trace(F'n)' — trace F'll' = — trace F'lT' = 2tracenF'F“^F', 

and 


2trace nF'^'«(F)F' - 2trace H F'F“^ ^s(F)F' - 2tracenF'F“^F' = -2tracenF'F“^“2^'. 
Hence (3.5) becomes 


d^Zjs) _ -r(s + l) 
diT^ r(s —n/ 2 ) 


2r(s +1) 


trace F"F“*“^ - trace nF'F“^“2 zp/ 


r(s — n/2) 


+ 


r(s — n/2) 


trace 


/o Jo 


{u + v)^F' - H) F' (e”"^ - H) dudv. 


We will reorganize this formula to combine terms which contain poles. We have 


trace 

trace 


u-\-v<l 


{u + v)^F' - H) F' (e“^^ - H) dudv = 


{u + v)^F' e-“-^ 


dudv — 2 trace 


(u + v)^F'U F'ie-^^ - m dudv 


fU+V<1 

= trace 


u-\-v<.l 




{u + v)^ F'F'dudv + 2 trace F'H F'F”! t 


dt - 


s + 1 


So defining L{s) as in (3.4), we obtain Lemma 3.1. 


4. Reducing to the heat kernel at small times. 

Our aim in this paper is to express the Hessian of Z{s) in the form given in (1.4). In this Section we 
reduce to a local problem involving the heat kernel at small times. Most of the steps in this reduction 
are quite routine. 
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We start by expressing F defined in (1.2) in local coordinates. Let be the cone of positive 

definite elements of There exist functions a, b, c, such that in the natural trivialization of E 

associated to arbitrary local coordinates on M, 

(4.1) F = + 6 + (9,5,c^/, 

ij i 

where I is the identity N x N matrix, and 

a : (0,oo) X 5+(M’") X r°’^(M’") ^ 
b : (0,oo) X Sl{R^) x 

c: (0,oo) X 52(M’^) ^ T^’^{R^)^T^'\R^) 

are smooth functions which are independent of the metric and the choice of local coordinates, so writing 
dg for the tensor {diQjk) G T°’^(M”) at each point of the coordinate chart, 

for 1 < i < n, a''= a^{V,g,dg) G T^’^{R^), 

h = b{V,g,dg) G T^’\R^), 
for l<i,j,kj<n, c^/ = c^/(y, 5 ) e T^’\R^). 

In the previous Section we computed the second variation of Z(s) when the varying metric is given 
by 

5 ( 0 -) = g + ah. 

Differentiating (4.1) once, we obtain 

(4.2) F' = ^ A% = A%{V,d'<g)^T^^\R^), 

l«l + i^l<2 

where 7 ranges over multiindices with \f3\ + I 7 I < 2. Differentiating (4.1) twice we obtain 

(4.3) 

F" = Yl = G T^’\R^). 

l«l + l/3|<2, 
l7l<2 

The Hessian of Z(s) with respect to the flat connection on the space of metrics is just given by 

(P 

HessZ(s)(/i,/i) = -^Zis). 

da^ 

We will show that if we can prove Theorem 1 in this case, then the result also follows when we use the 
intrinsic connection on the space of metrics. We also reduce (1.4) to the case when {{g, h))g = 0. We 
will then express the second variation given in Lemma 3.1, in terms of h, and analyze the terms which 
arise. 
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Definition. Write nj : M x M ^ M for the projection onto the factor. Fix a metric g on M. 
For a pseudodifferential operator Q : C°°{E) C°°{E) of order less than —n, define K{Q,x,y) to be 

the continuous section of tt^E 0 n^E* on M x M such that K{Q,x,y) dV{y) is the integral kernel of 
Q. If Q{s) is an analytic family of polyhomogeneous pseudodifferential operators then K{Q{s),x,y) 
may be analytically continued from the s-region where Q{s) has order less than —n, and we write 
K{Q{s),x,y) to denote this analytic continuation. When x ^ y, the section K{Q{s),x,y) is entire in 
s and K{Q{s),x,y)dV{y) gives the Schwartz kernel of Q{s) away from the diagonal. At the diagonal, 
K{Q{s),x,x) is a meromrphic in s. To compute the poles and the trivial zeros of K{Q{s),x,x), we 
refer the reader to [Okl] Sections 2 and 3. 

From (3.3) and (4.2), 

DZ{s){h) = traceF'L(s) = / trace K{F'L{s),x,x) dV{x) 

J M 

= f {d°^hij{x))K{A^^^d^L{s),x,x)dV{x). 

hj 

|a| + i/3|<2 

By introducing a partition of unity and integrating by parts, we can express this in the form 
(4.4) {{h,X{s)))g. 

where X(s) G The expression for L{s) in (3.4) was chosen precisely so that (3.3) holds and 

K{Ai^a<^L{s),x,x) 


is an entire function of s. 

Now we reduce Theorem 1 to the case when the following two conditions hold: 

(a) . We take the flat connection on A4 

(b) . {{h,g))g = 0. 

To make this reduction, assume that (1.4) holds when (a) and (b) are satisfied. We will first remove 
condition (b). For general h G C°°{S‘^{TM)), set 

g{a,T) = e^{g + ah), 


so 

5^g(0,0) = h, drg{0,0) = g. 

Then Z{g) is independent of r, so from (2.1) we see that 

Hessflat-2:(s)(/i,5r) = -DZ{s){d^drg) = - DZ{s){h). 

Define P : C^{S^M) C°°{S^M) by 

Ph = h- Chg, where Ch = ■\-{{h,g))g, 

nV 
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Then defining X(s) as in (4.4), 

Hessflat^(s)(/i,/i) = HesSflat-2i(s)(P/i,P/i) + 2c/, Hessflat 2:(s)(/i, 5) + Hessflat 2:(s)(c/, c/) 

= {{Ph,T,Ph))g - 2chDZ{s){h) 

= {{Ph,T,Ph))g - ^{{h,g))g{{h,Xis)))g 
= {{h,PT,Ph))g - {{h, ^{{g,h)) X{S)))^ 

The operator corresponding to the first term is PTgP which is a TDO equal to Tg up to the smoothing 
operator Tg — PTgP. The operator corresponding to the second term is /i ^ ■:^{{g,h))X{s) which is 
smoothing and analytic in s. This removes condition (b). 

By (2.2), for the intrinsic connection on A4 we have 

HesSintrinsic-2^(s)(^, ^) = Hessflat 2:(s)(/i,/i) + DZ{s){w) 

where 

Wij = ^ ^ hjkg hij + '^{h,g)g hij ~^{h,K)ggij 

k,i 

Then raising indices with respect to the metric g in the usual way, 

DZ{s){w) = - / V hik9^%gX^\s)dV + 11 {h,g)g{h,X{s))gdV - ] [ {h,h)g{g,X{s))gdV 

i,j,k,e 

The operator corresponding to this form is a degree zero operator which is analytic in s. 

This reduces Theorem 1 to the case when we take the fiat connection on A4 and assume that h is 
trace free relative to g. In particular we now have I^^(O) = 0. The rest of this section is devoted to 
analyzing the terms of the Hessian of Z{s) in Lemma 3.1, on the right hand side of (3.1). 

With L{s) defined in (3.4), from (4.3) we have 

(4.5) trace F"L(s) = f ^ hke){d^hij) trace L{s),x,x) dV{x). 

^ i,j,k,e 
|a| + l/9|<2, 
l7l<2 

The operator corresponding to this form has degree 2 and is entire in s. 

To deal with the remaining three terms in (3.1) we will need the following standard estimates on the 
asymptotics of the heat kernel K{e~^^,x,y). Write d for the geodesic distance between points x and 
y, and write V for the connection on E or E* associated to g. We take the norms | • \g associated to g 
on the fibers x M). There exists e > 0 such that for all j, /c, £, 

(4.6) sup d{x’^XyK{e~*^,x,y) < when t < 1, 

d{x,y)>do 

(4.7) sup (9/V^V^iL(e“*^ — n, X, y) < when f > 1, 

x,y^M 
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where the constant C depends on j,k,i. Moreover, if we write tt^ ; [0,1] x M x M ^ M, i = 1,2 
for the projection onto the ith M factor, and d{x, y) for the distance from x to y, there exists b G 
C°°{ttIE (g) 7r|£'*) and do > 0, e > 0 such that 


(4.8) K{e-^\x,y) = 


(47rt)”'/2i/ 




when t < 1, d < do- 


Moreover, b{0,x,x) is the identity matrix in E*. See for example [Ch], 

By (4.7), we see that 

/ OO 

— n, X, y) dt 

is a smooth kernel which is entire in s, and so there exists a smoothing operator P{s) which is entire in 
s such that 


(4.9) trace (^F'UE'F-^ P {e-^^ -U)dt^ = {{h,P{s)h)). 

Now consider the term from (3.1): 


dudv 


(4.10) 


trace I I {u + vf F' (e““^ - U) F' (e“"^ - U) dudv 

/ u-\-v>l 

n oo 

” 

^1/2 ^oo 

+ trace / / ” 

J u=0 Jv=l—u 

/»l/2 poo 

+ trace / / ” 

Jv=0 Ju=l—v 

= (I) + (II) + (III). 


dudv 

dudv 


In any of these terms we can pass the trace inside the integral. Indeed, writing || • || for the operator 
norm, || • ||i for the trace class norm and and || • II 2 for the Hilbert-Schmidt norm of an operator, when 
u > 1 / 2 , 

||F' (e-“^ - n) F' {e-^^ - n)||i < ||F' - n)|| 2 ||(e -“^/2 - B) F'|| 2 || (e”^-^ - n)|| 

< Ce-™, 

where the last inequality follows from (4.7). This shows for example that 

^ 1/2 poo 

trace / / (u + z;)®||F' (e““'^ — H) F' — n)||i dudr; < 00 . 

Jv=0 Ju=l—v 

so that the trace can be passed through the integral in (III). 

Now we localize by writing 1 = YlxeA ^ where A is a finite partition of unity on M such that for 
each point x G M, all the functions x G A which contain x in their supports are supported in a single 
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contractible coordinate chart, over which E is trivial, and on which d°‘gij, and are uniformly 

bounded for each a. Hence we express 

(4.11) trace F' - H) F' - H) 

= ^ trace - H) xA’^^ {d^hki)d^ - H) . 

XiXGA, i,j,k,£ 

|a| + l/3|<2, |7l + |5|<2 

We can analyze each term in the above sum separately. Fix XTX:hj:k,£,a,f3,x,6 and set 

(4.12) E{u,v,x,y) = x{x)x{y) trace A^^^{x) {d^K{e~^^ -U,x,y)) A^^^iv) {dyK{e~^^ - U,y, x)). 
The trace on the right is the trace on N x N matrices. Then 

(4.13) trace (e'^^ - H) {d^hM)d^ - H) 


/Supp(x)xSupp(x) 

Now from (4.7), for tt > 1 and u > 1 and each a, f5, 

\d^d^E{u,v,x,y)\ < 

Hence 

1 


{d^hij{x))E{u, V, X, y){d^hki{y)) g{x)^g{y) dxdy. 


(4.14) 
is smc 

(4.15) 


KP{x,y) := 


noo /*oo 


{u + vyE{u, V, X, y) dudv 


r(s-n/2 )7i 

is smooth in (s, x, y) and entire in s and (I) is a finite sum of terms of the form 

{d°hij{x))KP{x,y){d'^hkeiy)) ^g{x)^g{y) dxdy. 

Supp(A) xSupp(S) 

Integrating by parts to take the derivatives from h, we find that the kernel of this bilinear form is 
smooth and entire in s. Similarly, write (H) formally as a sum of terms of the form 

(4.16) ff (d^^hijix)) (^x{x) traceiFp)(x,y) xiv) ) hkt{y)) \/g{x)y/g(y) dxdy, 

J t/ X O 

where Ka^^\x,y) is given by 

/» 1/2 poo 

Ki^^\x,y) = / / (u + u)M*J^(x)(afF:(e““^,x,y)) yl^l(y) (a^iF(e“^'^,y,x))dudu. 

Ju—0 Jv=l—u 

In doing this we are switching the dxdy integration with the dudv integration, which must properly 
interpreted and justified, since the integral in (4.16) may be divergent due to Ks^^\x,y) being singular 
on the diagonal. The way to proceed is to work in the class of pseudodifferential operators of order \f3\. 
Define 


4^^\u,x,y) := A^liy) 


f V=-l — U 


u + vY{dlK{e ,y,x))dv A"^^{x) 


so 


(4.17) 


fl/2 


Ki^^\x,y) = / {d^K{e ,x,y)) jy^\u,x,y) du. 


/m=0 


From (4.7), Js^^\u,x,y) is smooth for {s,u,x,y) G C x [0, |] x D x D and analytic in s. We notice 
that the integrand on the right hand side of (4.17) is the integral kernel of an operator on C'“(D) 
whose order \f3\ symbol norms are bounded uniformly in the parameter u. This follows easily from the 
elementary lemmas below, and shows that the integral in (4.17) converges to the Schwartz kernel of an 
order |/3| pseudodifferential operator. 
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Lemma 4.1. The operators e are psedodifferential operators of order zero, with symbol norms 
uniformly bounded in the parameter u. 

This implies that for j3 fixed, the operators are psedodifferential operators of order |/3| with 

uniformly bounded symbol norms in the parameter u. 

Lemma 4.2. If K,{Q,x,y) is the Schwartz kernel of a 'ItDO, Q, on M of order m, and J{x,y) is 
smooth. Then the product K,{Q,x,y)J{x,y) is the Schwartz kernel of a It DO, which we denote by Qj, 
which also has order m. If Q depends smoothly on a parameter p £ then so does Qj. 

We can conclude that the expression in (II) is equal to (4.16), and Ki^^\x,y) is the Schwartz kernel 
of a TDO of order \f3\. The question then arises whether this operator is polyhomogeneous and how 
one obtains its symbol expansion. For a function ■0 £ ^“([0,1]), integration by parts gives 


fl/2 


e-^^f^{u)du = A-V(O) + A-20'(O) + ... + 


/ti=0 


- (a-V( 1) + A-20'(l) + ... + A-^0(^-^)(l)) e-^/2 + J ^ e-^^f^^^\u)du, 

we see that (4.17) equals 
(4.18) 

M-l M-1 

^ (5fiF(F“’^“0x,y)) dff\u=oJi^^\u,x,y) - ^ (^f x, y)) dff\u=iJi^^\u,x,y) 


m=0 


m=0 


. 1/2 


+ 


(af F(F-^e-“^, X, y)) 4") (u, x, y) du. 


lu=0 


The term in the second sum 


-{d^K{F ^ ^e ^/^x,y)) 50"|„=ijp)(u,x,y) 

is smooth in (s, x, y) and analytic in s. We need the following refinement of Lemma 4.2; 

Lemma 4.3. In fact if the It DO, Q, in Lemma f.2 is polyhomogeneous of degree fi, then so is Qj, and 
if ,1 = Ja depends analytically on a parameter s, then Qj^ is an analytic family of pseudodifferential 
operators. 

With this refinement, we see that the term in (4.18) 

(9f iF(F— 1, X, y)) {u, X, y) dV{y) 

is the Schwartz kernel of a polyhomogeneous ikDO operator of degree \f3\ — 2m — 2. The last term in 
(4.18) is the Schwartz kernel of a pseudodifferential operator of order \f3\ — 2M — 2. We conclude that 
Fi") (x,y) is the Schwartz kernel of a polyhomogeneous ikDO of order |/3| — 2, which is analytic in s. 
The lemmas are consequences of the standard theory. Indeed for Lemma 4.1, following [Gi] one can 
express the symbol of the heat operator in terms of the symbol of the resolvant via a contour integral. 
However, the resolvant (A — F)~^ is a pseudodifferential operator, and the zeroth order symbol norms 
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of A(A — F)~^ are uniformly bounded in A. Lemmas 4.2 and 4.3 follow from the relationship between 
amplitudes and symbols, see [Tr]. 

We return to our analysis of (4.10). Term (III) is dealt in the same fashion as Term (II) with the 
roles of u and v interchanged. Now we want to understand the final term in (3.1) which is given in 
( 1 . 10 ): 

(4.19) ——^^ trace f f (u + vYF' F' e~'"^ dudv. 

T{s-n/2) JJu+v<i 

From now on we will refer to the steps listed in Table 1.1. 

Step 14: Using estimates based on (4.8), we can take the trace in (4.8) inside the integral because 
||F'e““-^ F'< ||F'F'||2||e“^-^|| < 

Since this is also true with u and v interchanged, we have 

||F'6-“-^ F'< C{u + v)-^-\ 


and so for sufficiently large, 


(u + u)'*||F'e “ F'e " \\idudv < oo. 


u-\-v<. 1 


Steps 7, 3: Taking a partition of unity 1 = J^xeA X on M as for (4.11), we see that 


(4.20) trace F'F' 


^ trace y A% id^hij)d^ 6““^ X {d^hu)d^ e 


-vF 


X,x€A, 

|a| + l/3|<2, |7l + l<5|<2 


Step 1: Fixing x,x G A, i, j,/c, £, a,/3,7 , d, we will consider each term in (4.20). Set 

(4.21) E{u,v,x,y) = x(x)x(y) trace A(^^^(x) (^fF(e““-^, x, y)) A^|( 2 /) {d^K{e~^^,y,x)). 

Step 10: Then 

(4.22) trace {d^hki)d^ 


'Supp(x)xSupp(x) 

Steps 8 , 14: Hence formally, (4.19) is a sum of terms of the form 


(9“%(x))F(u, V, X, y){d'^hkeiy)) g{x)^/g{y) dxdy. 


(4.23) 

where 

(4.24) 


Supp(x)xSupp(x) 


Ks{x,y) := 


{d°'hij{x))Ks{x,y){d^hM{y)) \/g{x)g{y) dxdy, 


r(s - n/2) JJu+v<i 


{u + vYE{u, V, X, y) dudv. 
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When the supports of x and x are disjoint we can use (4.6) and (4.8) to see that (4.19) is a bilinear 
form on {S‘^M) with a smooth kernel which is analytic in s. Hence we have reduced to the case 
when X and x are supported on a coordinate chart H whose diameter is less than the injectivity radius 
of M. 


5. The symbol of Ug- 

Here in Section 5, we describe how to compute the symbol expansion of Ug, and we carry this out 
explicitly for Ug thus obtaining Theorem 2. To do this we expand on the summary in Section 1. Section 
5 can be read independently from Section 6, which contains the rigorous proof of Theorem 1, and shows 
that the algorithm we lay down here is correct. We will see in particular in Section 6, that (4.19) is a 
sum of terms of the form (4.23). Moreover, we see in Lemmas 6.1 and 6.5 that Kg defined in (4.24) has 
symbol 

Ugix,0 + 14(x,0, 


where setting 


m = \p\ + |5|, 


Ug is poly homogeneous of degree n — 2 s — 4 + m, and Vg polyhomogneous of degree m — 2 . 

Step 3 : We are interested in obtaining the symbol expansion of Ug at the point x G M, and so we 
work in coordinates on a small coordinate chart H containing x. 

Steps 11, 12, 13: To obtain the symbol expansion of Ug, we make an expansion of E{u, v, x, y). Using 
(4.8), 


(5.1) 


d^K{e-^^ ,x,y) = 


where x, y) is a smooth N x N valued function on [0,1] x whose Taylor expansion at u = 0 

gives an expansion of the form 


b(i{u,x,y) ~ ^ uH 0 i{x,y), bfsi{x,y) = 0(d(x,y)”^^^ l/3|,o})^ 

i=-m 

Here, b 0 i{x,y) is smooth on Hence considering (4.21), 

(5.2) E{u,v,x,y) = J{u,v,x,y), 

where 

J(^, y-^ 

is smooth on [0,1]^ x and J{u,v,x,y) has an expansion of the form 

(5.3) 

oo oo 

J{u,v,x,y) ~ E E u‘'y'Jab{x,y), Jab{x,y) = 0(d(x,y)“""^-2(“+^)-(l^l+l^l)’°>). 

a=-|/3| 6=-|c5| 

Here, Jab{x,y) is smooth on 
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Step 4: Fix X and write w for geodesic normal coordinates of y centered at x . We make a Taylor 
expansion of Jab{x,y) about re = 0 to write 

OO 

(5.4) J{u,v,x,y) ~ 'h{u,v,x,y), 

i— — m 

where for constants Jabnix), 

(5.5) Ji{u,v,x,y) = ^ Jab^lU°‘v'w^. 

2 a+ 2 b+\^l\=^ 


The sum in (5.5) is clearly finite. Now set 

(5.6) 

K.A^;y) = r(, - „/ 2 ) 

Step 15: The symbol of E£(u,v,x,^) has the form 


Ei{u,V,X,y) = 

1 /■/■ 

(u + vyEi{u, V, X, y) dudv. 


Us,yx,i) + Vs,i{x,y), 

where Us/{x,^) is homogeneous of degree n — 2s — 4 — and Vs{x,^) is polyhomogeneous of degree 
m — 2. To see this, the metric-modified symbol (see 2.4) of Ei{u,v, x,y) is obtained by taking the 
Fourier transform in y, 

(5.7) a{Ei,u,v,x,^) = J e™'^Ee{u,v,x,y) dw. 

From the definition of E^, 

a{E,,u,v,x,0 = (47r)"/2(n + n)-"/2 

2 a+ 2 b+\bi\=^ 

(5.8) 

2 a+ 2 fo+|^|=£ 


= (n + n)-^/2-^g-|?lW(«+N J„b^(x)u“nV, 

\v'\<2m-\-£ 

2a-\-2b=3i-\-\iy\ 

where the constants c^i, and functions Jabv{x) are defined appropriately by collecting terms. We will 
prove in Lemma 6.3 that Jabw{x) vanishes if a or 6 is negative, so 

a{Ei,u,v,x,^) = 0((n + n)“"'/^“^) 


24 



as u,v ^ 0 . 

Step 14: Hence the symbol of Ks/{x,^) in (4.24) equals 


(5.9) 


1 

r(s — n/ 2 ) 


JabAx)C + 

\u\<e JJu+v<i 

2a+2b=3£+\u\ 


which converges absolutely for 3?s large. 

Step 16: By making the change of variables 

T=iu + vM\\ 


T = 


U 

u + v' 


we see in Lemmas 6.4 and 6.5 (also (1.19)) that this integral extends analytically and to a sum of 
symbols of degrees m — 2 and n — 2s — 4 — i. The symbol of degree n — 2s — 4 — i is obtained by 
extending the integral on the right of (5.9) to {u,v) G [0,oo)^, and so 


(5.10) Us,eix,0 = E JabAx) C 

w\<i 

2a+2b=3e+\u\ 

r(s — n/2 — i + a + b + 2)r(n/2 — s + i — b — l)r(n/2 — s + i — a — 1) 
r(s — n/2)r(n — 2s + 2.^ — a — 6 — 2) 

One might hope to complete the analysis of term (4.19) by studying the properties of the remainder 
term in the Taylor expansion (5.4). However although the order of Us,i decreases as i increases, the 
order of Vg/ does not. We will need to decompose the kernel E{u, v, x, y) in additional ways which pick 
out the homogeneous terms in the symbol of Vg- This is carried out in Section 6. We emphasize that 

(5.10) gives the term of degree n — 2s — 4 — £ in the symbol of Kg{x,y). From (4.23), Us{xA) is the 
symbol of 

(5-11) E (-1)'“'^''^' > 

X.xeA, i,j,k,£ 

|a| + |/3|<2, |7l + l'5|<2 


which can be computed from (5.10) using the symbolic calculus. It is also easy to go from the metric- 
modified symbol to the true symbol using the symbolic calculus if one so desires. See [Tr] for details 
on the symbolic calculus. One needs to be a little careful because we have been working in normal 
coordinates based at x, and in order to apply the usual symbolic calculus formulas to compose operators, 
one has use the same coordinate system for the points x and y. The necessary change of variables 
formulas can also be found in [Tr]. 


Proof of Theorem 2. We want to compute Ug{x,f/) which is the degree n — 2s term in the symbol of 

(5.11) . (As always, in saying this we must assume s is not a half-integer greater than n /2 — 1 , since at 
these values the symbol of H(s) also contributes to the degree-(n — 2s) term.) Now Ug^-m{x-,f) is the 
degree n — 2s — 4 + m term in the symbol of Kg{x, y), which is the principal symbol when is small. 
Hence from (5.11), 

(5.12) n,(x,0 = E (-l)l“l(i6“+^t^«,-m(x,0- 

|a| + |/3|<2, |7l + F|<2 

We will evaluate 5.12 at a fixed point x = xq, and from now on we express x and y in geodesic normal 
coordinates centered at xq. In particular, we assume that the derivatives in (4.21) are computed in 
these coordinates. 
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Lemma 5.1. Let C{s) be defined as in (1.7). Define 


(5.13) G = ^-m^uv/{u+v)^ 

define the polynomial P{d°',t,w) in the variable w G M"' by 

(5.14) P{d^,t,w) = 
define 

(5.15) cj(5“, d^, d\ u,v,x,() = (iO“+^ P{d^,u,-id^) P{d^,v,-id^)G, 
and define 

1 


(5.16) Us{d^, 5^, dfi x,C) = 
Then 

(5.17) Us{xo,fi) = G(s) 


n oo 

(u + t')®(j(5“, , 5'*', d^, u,v,x,fi) dudv. 


G{s)T{s-n/2) 7o 


^ «s(<9“, d^, d'^, d^, xo,Oti'aceyl*^^^(xo)^!^5(xo). 


|a| + |/3|=2, |7l + |5|=2 


Proof of Lemma 5.1. Recall that each term Ks{x,y) in (4.24) depends on the parameters XiX ^ 
A, i, j, k, a, /3, 7 , <5, and is given by 


Ks{x,y) := 


1 


r(s — n/ 2 ) 


(u + vyE{u, V, x, y) dudv. 


where 


E{u,v,x,y) = xix)x{y) trace A"^^{x) {d^K{e , x,y)) A')l{y) {d^yK{e ^^,y,x)). 

The degree n — 2s — 4 + m term in the symbol of Kg is Ug^-mix,^), where 

pOO poo 

(5.18) Ug-m{x,fi) = ^ / / {u + Vy(T{E_rn,U,V,X,f)dudv, 

r(s - n/2) Jo Jo 

and E-rn which gives the lowest order terms of E is defined in (5.6). We will see that writing w for the 
geodesic normal coordinates of y about rc = xq, we have 

(5.19) E_^{u,v,xo,y) 

= xixoMxo) ( 47 r)-/ 2 (uu)-/ 2 (-l)l/ 5 |(a^e-l-l^/ 4 “) traceA^^(xo)A^i(xo). 

To see this, from (5.1) we must compute the “lowest order” terms in the functions 

P (xo ,y ) 14:V r.6 

, C L/y 

X=Xq 


e-"'("-")/4”6(u,xo,y)) 


^P{x,y)l^uQf3 (^e-'^'("’^)/4“6(n, X, y)) 
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about u = V = 0 and y = xq, where the orders of the terms u°"w^ and v°'w^ are 2a + \y\. It is clear 
that to get the lowest order terms, we can replace b{u,x,y) and b{v,xo,y) by 6 ( 0 , xo,xo) which is the 
identity matrix. This reduces the problem to computing the lowest order terms in the expansions of 
the functions 


(5.20) 

(5.21) 


^d'^{x,y)/4uQy^-d'^(x,y)/4:U\ ^ 


X=Xq ’ 


^d'^{xQ,y)lAvQ5^-d'^{xo,y)l4v ^ 


Now since the derivatives are expressed in geodesic coordinates about xq, (5.21) is precisely 

\w\^lAv nS -\w\^jAv 
c. 1 


which is a linear combination of terms v°‘w^ with order 2a + |;u| = —1(5|. On the other hand (5.20) may 
involve terms coming from the curvature of the metric. However, these are all of higher order than 
—1/3|, as one can check by expressing x and y in normal coordinates around xq, and expanding 

d{x,yf = |x-yp + O {\x - y\^) . 

We thus find that the lowest order terms in (5.20) are given by 

g|a:-yP/4u QP^-\x-y\^/Au I _ qI 3 / Au 

From this, we get (5.19). Since the Fourier transform of is G, taking 

the Fourier transform of if-m gives 

(5.22) a{E_^,u,v,XQ,i) = y:(xo)y(xo) traceH*J^(xo)H^ 5 (xo) (-1 )I^Ih(5^,u,- i5^) P(5'^,x, 


Plugging this into (5.12) and (5.18) and summing over x,x G A yields Lemma 5.1. 

To complete the proof of Theorem 2, it remains to compute u(cl“, 9^, 9'^, 5'*', u,v,x,^) for explicit 
values of a, f3, 7 , 5. Now 


Wi 


P{dj,u,w) = 


2u' 


p{dj,u,-id^) = —d^^. 


P{djdk,u,w) = 


djk 

4^2 2u 


P{djdk,u,-id^) = 


= -(P 


\ 2 %% + ^ 

l 4^2 ‘ik 2 u 


and 


a G = 

^ u + v 

4 u^x^ 


= 


{u + v) 


idik - ] G, 


2^1 


U + V 


8 u^v" 


■,3 Au'^v^ 

= [ ~777T++3^3^k^p + I „,^2 (djk^p + bjp^k + bkpij) ] G. 


{u + u)" 


8 u^v^ 


{u + v)^ 
4 u^v‘^ 


{Sjkbpq T bjpSkq T bkpbjq) 


{u + v)^ 


{u + x)^ 

(.djqCkCp T dkqCj^p T dpq^j^k + bjkCpCq T ^jp^k^q T ^kp^j^q) ] G. 
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So for (a), 
for (b), 


cri^djdk, /, dpdq, /, Cj^k^p^qG 


a{dj, dk, dpdg, I, u,v,x,C) = -iCjCpCqPidk,u,-id^)G = =- ^ijikipiqG. 

ZiU U Z) 


for (c), 




a{dj, dk, dp, dq, u,v,x,^) = -^j^pP{dk,u,-id^)P{dq,v,-id^)G = 


2^j^kCp^q I ^^^^kqCjCp ] G. 


{u + v) 


for (d), 


_ f 


a{I, djdk, dpdq, I, u,v,x,C) = -^pCqP{djdk,u,-id^)G = + 


2{u + v) 


^p^q^jk A 
2u 


G 


4u^ \{u + v)^ u + v 

for (e), 


2u 


G = - - ^ijikip^qG + 7r7~~~:^^jkCpCqG- 


CpCq^jk « 'V 


{u + v) 


2{u + v) 


a{I, djdk, dp, dq, u,v,x,i) = iipP{djdk,u,-id^)P{dq,v,-id^)G 

= fJ^d d.,d. + ^d. IG 

\8u^v Auv 


?/?y ?) 1 

-^ijikipiqG + r,^„. I -.02 i^jk(p(q + ^jqCk^p + ^kqCjCp)G - ^ __^ SjkCp^qG 


(u + v) 


2{u + v)^ 


2{u + v) 


and finally for (f), 

a{I, djdk, I, dpdq, u,v,x,Cj = P{djdk,u,id^)P{dpdq,v,id^)G 


u^v'^ 


'CjCkCp^q + 


1 


{u + vY ^ p 9 ' _i_ 

uv 


{5jkSpq p 5jp5kq P ^kp^jq) 


2(u P v)^ 


i^jqCkCp “1“ ^kqCjCp P ^pqCjCk P ^jkCpCq + ^jpCkCq + ^kpCjCq) ) G 


U \ V 11 

+ ( 7:7 ^ —^jk^pCq — Tr~T~P~.^jk^pq + r,f„, , ~ A„.(.. i „.\^i^^Pq ^ ^ 


2{uPvY ^ P 1 4:{uPv)v' 

1 


2{u + v) 


Au{u P v) 


u^v'^ 


p 


[u P v) 

V 


^ijik^p^q + 


2{u P V 


2 ij^kSpq 


A{u P vY 
uv 


{5jkSpq p Sjpdkq P 5kp5jq) + 


u 


i^jk^p^q 


2(u P v)^ 


2(n + ^;)2-J-'’P'’9 

ipjq^k^p + ^kqCjCp P ^pq^j^k P ^jk^pCq P ^jp^k^q + ^kp^j^q) ) G. 
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Now we just need to substitute these results for a into (5.16) to obtain the expressions for Us given in 
Theorem 2 (a)-(f). This is done more efficiently by looking at each type of term which occurs in a and 
plugging it into (5.16). Writing S = s — n/2, we have 


1 


r^OO /‘OO 


r(s-n/2)yo Jo 


{u + vf-^uPy^^Gdudv = 


/ 2 .^- 2 S+ 2 r- 2 p- 2 <,- 4 '^(.S - r + p + q + 2)r{-S + r - p - l)r{-S + r - q - 1) 
^ ^ r{S)r{-2S + 2r-p-q-2) 


S---{S+l+p + q-r) {2S-1)---{2S+ 2+p + q-2r) 
S ...{S+ l+p-r) S---{S+l+q-r) 


cm\ 


-2S+2r-2p-2q-A 


where C{s) defined in (1.7), and we assume r < min{p, g} + 2. Hence we obtain Us{d°, d^, , 

from (t(9“, 5^, 5*^, 5'*', u,v,x,^) by substituting 


^ ( 5 ^ + 5 ) |er---^ 

^ ( 25 ^ + 5 -l)|Cr^^-^ 
{AS^ + 2S - 2) 

[u + v)^ 

^^G ^ (45^-25)ler---^ 

G 4(4,5^ - l)|^p-2"-4. 


uv 


{u + v)^ 

V 

{u + u)2 

1 


G 

G 


-G 


U + V 

1 

{u + v)‘ 


:G 


(25-i)ier“"^ 

(4,S-2)|^|^-2"-2 

- ler- 2 ^ 


6. The double decomposition. 

In this Section we complete the proof of Theorem 1. In order to simplify the notation we will restrict 
to the scalar case with V = 1. It is straight forward to obtain the vector case by replacing scalars 
by matrices in the appropriate places. Let 17 be a coordinate chart of M of small diameter which is 
identihed with a subset of M"", and make a change of coordinates {x,y) {x,w), sending Q x Q onto 

A C where {x,y) are product coordinates on 17 x f7 and w = w{x,y) are normal coordinates of y 
centered at x, which vary smoothly with x. We define the function p(x, w) to give the volume element 
on M; 

p{x,w)dw = dV{y)= y^g{y) dy. 

For a function / on 17 x f7, we define /•^ on A by 

f{x,y) = f^{x,w{x,y)). 

and we are interested in understanding (4.19), which is a sum of terms of the form 
r(s - n/2) JJu+v<iJ Jnxn 
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where 


(6.2) E{u,v,x,y) = A{x){d^K{e ,x,y)) B{y) {d^K{e '’^,y,x)). 

Here, the smooth functions A and B are supported in H. Set 

(6.3) m = \/3\ + |(5|. 

Step 9: We denote the operator on H with kernel E{u, v, x, y) by £u,v and the metric-modified symbol 
of this operator (see (2.4)) by a{E,u,v,x,^), so 


(6.4) a{E,u,v,x,^) = f e™'^E-^{u,v,x,w) dw, 

and for / G C°°(Q) compactly supported, 

f cr{E, u, V, X, ^)e“™'^/'^(x, w)p{x, w) dwd^. 


{£u,vf)ix) = / E{u,v,x,y)f{y)^/g{y) dy = 


(2vr) 


Steps 14, 15; We will see in Lemma 6.3, that 

\a{E,u,v,x,0\ = + 


as u, u —> 0. 


Hence 

i£u,vf){x) = 0((u-h asu,i;^0, 

and when > n/2 -|- m — 2, we can interchange the {u,v) integration in (6.1) with the x integration 
and conclude that (6.1) equals 


(6.5) 


where Ws is the operator 


which has symbol 


[ {d‘^hij) {Wsd^hki) y/gdx. 
Jn 


T{s-n/2) JJu+v<i 


{u + v)^£u,v dudv, 


( 6 . 6 ) 


crs(x,0 := 


r(s-n/2) JJu+v<i 


{u + vYa{E, u, V, X, ^) dudv. 


Our aim is to show the following Lemma. 

Lemma 6.1. as(x,^) is entire in s and has the form Ug -|- 14 where Ug has degree n — 2s — A + m, 14 
has degree —2, and Ug and 14 axe meromorphic with simple poles in n/2 -|- N"*". 

This completes the proof of Theorem 1, since we can now write (6.5) in the form 


hij 


i'Tghki) ^/gdx. 
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where 7^ is the sum of poly homogeneous ^'DOs of degrees n — 2s and 4 — m — 2 < 2. We pause to 
outline the strategy to prove Lemma 6.1. It is convenient to make the change of variables 

(6.7) t = u + v, T = u/{u + v), 


so 

Write 


u = tr, 


V = t{l — t). 


( 6 . 8 ) 


H{t,T,x,y) = E{tT,t{l -T),x,y), 


so 

(6.9) 


a. 


= 


r(s — n/2) 


t, T, X, drdt. 


0 JO 


By (4.8), we see that 

(6.10) H{t,T,X,y) = t-^-"^(r(l-T))-”/2-mg-d^(x,y)/(4tr(l-r))p(^^^^^^y) 
where 

(6.11) P{t,T,x,y) G ^“([ 0 ,1]2 X 

Steps 12, 13: We make a Taylor expansion of P{t,T,x,y) = P-^{t,T,x,w) in w. The order 
remainder term has order \w\^, and in Lemma 6.2 we see that it contributes a symbol of order —L + C 
to iTs(x,^), where C = C{m,n,s) is independent of L. The consequence of this is that we can ignore 
the Taylor remainder term and prove the result when P{u, v, x, y) is replaced by a Taylor polynomial. 
Step 16: Doing this and carrying out the r integral in (6.9), we find in Lemma 6.4 that cJs(x,^) is a 
finite sum of terms of the form 

( 6 . 12 ) C [ P^^Qit,\C\H,x)dt, 

Jo 

where Q{t, T, x) is smooth on [0,1] x [0, oo] x D (including the point T = +oo), and vanishes at T = +oo. 
We show in Lemma 6.5, that such a symbol is entire in s. Making a Talyor expansion around t = 0 and 
T = oo enables us to obtain the polyhomogeneous expansions of Ug and 14. Terms of increasing power 
in t give terms of decreasing order for Ug, while terms of decreasing power in T give terms of decreasing 
order for 14. We proceed to carry out this plan. 

Proof of Lemma 6.1. 

Lemma 6.2. Suppose that Og is given by (6.9) where H is given by (6.10) and P satisfies (6.11). Fix 
s gC, and /c G N, and suppose 

P'^ {t,T,x,w) = 0{\w\^) 

where 

L > k + {n — 2iR.s — 2)_|_ + 2m — 2. 
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Then as{x,^) is a symbol of order —k. 

This lemma is proved later on. It is a standard fact in the theory of 'I'DOs that given a a sequence 
of symbols on Q with Sj{x,f) homogeneous of degree fi — j, then there exists a symbol 
with symbol expansion 

OO 

i=i 

Moreover, if each term cjj depends analytically on a parameter s, then a can be chosen to depend 
analytically on s. This combined with Lemma 6.2 enables us to restrict attention to the case when 
P-^{u,v,x,w) is replaced by its Taylor polynomial in w of degree L; 

PLit,T,x,w) := ^ Pf,{t,T,x)w^". 

\n\<L 


Steps 12, 13: Set 

(6.13) HL{t,T,X,w) = (r(l - r))-"/2-mg-|«;|V(4tr(l-r)) ^ P^{t, T, x)w^^. 

\n\<L 


Step 15: Then 

a{HL,t,T,x,(,) = [ e^'^'^HL{t,T,x,w)dw 

Jkp 

= t-^-"^(r(l-T))-’^/2-m ^ P^{t,T,x) [ du; 

M<L 

M<L 

(6.14) = Pf,{t,T,x) ^ (tr (1 - t) ) V, 


M<L 


|/i.+i/| even 


where c^jy are constants. The problem with this expression is the presence of the non-integrable factor 

(r(l-T)) —, 


which may not be neutralized by the factor 

(t(1 -r))(l'^+'^l)/2 

when \n\ is small. However, the following lemma shows that this is not an issue after the terms are 
summed up. We see from (6.14) that we can write a{P[L,t,T,x,^) in the form 

a{HL,t,T,x,0 = ^ r, x)^• 

W\<L 
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Lemma 6.3. If L > 3m, then 


G C°°([0,1] X [0,1] X Q). 


This lemma is proved later on. We now need to improve our understanding of the dependence of 
a{HL,t,T,x,f,) on t. Indeed, Paying closer attention to what happens when we differentiate the heat 
kernel, we see that (6.13) can be improved to 


and so from (6.14), 

Hence writing 
we see that 

RL,u{t,T,x) := 
Now up to a term of order —L + C, 


P^it,T,x) = 0(t("^-|^l) + /2), 

RLAt,r,x) = 

[r] := inf{g G Z : q > r}, 

G ^“([O,!]^ X H). 


(6.15) <Ts(x,0 = 


^ V ' 



(.-n/2+r(kl—.)/2l + l,,-|S|’>r(l-T)j;^^^(,_ T, x) 


0 Jo 


It is convenient to compactify the interval [0, oo) by including the point +oo. The differentiable structure 
at +00 is given by the coordinates f : (0, +oo] ^ [0, +oo). 




l/T T < +00 
0 T = +00 


Lemma 6.4. If 


R{t,T,x) G C°°([0,1] X [0,1] X Q). 


then ^ 

Q{t,T,x) := j e~'^'^^^~'^A{t,T, x) dr 
Jo 

is smooth on [0,1] x [0, +oo] x H and vanishes atT = +oo. 

This lemma is proved later on. We need one more result to complete the proof of Lemma 6.1. 
Lemma 6.5. Suppose Q{t,T,x) is smooth on [0,1] x [0,+oo] x H and vanishes atT = +oo. Then 

a{s,x,f) := f P~'^^‘^Q{t,\f,\‘^t,x)dt 
Jo 

is a symbol of order —2 on the set iR-s > n/2 — 1, which extends analytically to a meromorphic function 
for s G C, and 


33 



is smooth in {s,x,^) and entire in s. Moreover, 


cr(s,x,0 = Us{x,^) + Vs{x,C), 


where 


1 

sin(7rs) 


Us{x,0, 


1 

sin(7rs) 




are poly homogeneous of degrees n — 2s — 2 and —2 respectively, smooth in (s,x,^) and entire in s. The 
symbol expansions of Us and Vs can be computed in terms of the Taylor series of Q. Indeed, if 


Q{t,T,x) ~ '^t’"qk{T,x), t^O, 

k=0 

oo 

Q{t,T,x) ~ T^oo 

e=o 


then as |.^| ^ oo, 


Us{x,f) 

ysix,o 


oo « 0 O 

Y,^s,k{xMr-^^-^-^\ where Us,k{x) = / T^-^/^+>^qk{T,x) dT, 

k=0 

oo „1 

'^Vs,k{xM\~‘^~'^'", where Vs,k{x) = / U~'^/‘^~^~^q^{t,x)dt. 


We apply Lemmas 6.4 and 6.5 to (6.15). In order to apply Lemma 6.5 we replace s by s + — 

m)/2] +1 and deduce that (Ts(x, is analytic for s > n/2 + m/2 — 2, and is a sum of polyhomogeneous 
symbols of degrees —2 and n — 2s + m — A which are meromorphic in s with poles in n/2 + N"*". Since 
m is an integer between 0 and 4, this completes the proof of Lemma 6.1. □ 

The rest of this paper is devoted to proving Lemmas 6.2, 6.3, 6.4 and 6.5. 

Proof of Lemma 6.2. We consider the kernel Ks with symbol 0 "^. This is given by 

Ks{x,w) := ——[ [ U^^Hy{t,T,x,w)dw 
T{s - n/2) Jo Jo 

= ^ , , [ [ P{t,T,x,w) drdt. 

T{s - n/2) Jo Jo 

Now suppose K is a multiindex and |k| = k. Then 


\d:sd:/Ks{x,w)\ 
1 


1 rl 


T{s- n/2) Jo Jo 
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Using Leibnitz rule for differentiation, we see that this is bounded by a linear combination of terms of 
the following form where (. is in the range ^ < ^ < k. 


1 fi 


w 


L+2l-k 




0 JO 


and setting T = \ w\ this has the form 


\w 


r\w\ ^ 

L — k-\-2'?i^s — 2n — 2m-\-4: I — n — m — i I 


(r(l - drdT. 


>0 


>0 


Now 




T ^ 0 
T ^ oo, 


so we conclude that 

\d^d^K,{x,w)\ < C(n,s)|u;|^-'=+2K.-2^-2”*+4(l + |u>|”-2^"-2log(c/|u;|)), 
where c is a constant larger than the diameter of M. This implies that d^d'^Ks{x, w) is integrable when 
(6.16) L > k + {n-2^s-2)+ + 2m-2. 

Now suppose that (6.16) holds for some integer k > 0. Then 

as{x,^) = f e''^'^Ks{x,w)dw, 

and 

d^d^as{x,^) = f e^'^'^{iwYd"Ks{x,w)dw. 

JKX 

For fixed there exists a j such that 


101 > 


10 


n 


Then define the multiindex k, by 


= 0 £ / j, 

Kj = k + |/U|. 


so 


in > where C = ^(n, A:, ^u) > 0. 


Since 


e™-^ = 


(^ 0 ' 


-(9Ce™-€), 


we can integrate by parts to obtain 

d^d^asixY) = , f d'^d^{{iwYKYx,w)) drdtdw. 

V *sj JR’» 
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Replacing L and k in (6.16) by L + \ij,\ and /c + |/i|, and choosing R greater than the diameter of M, 


The integral on the right is finite when (6.16) holds. This completes the proof of Lemma 6.2. □ 

Proof of Lemma 6.3. We will start by listing the relevant definitions. Set 

A = {(x, re) : (x, y) G C 11 x R"". 

We have 

E{u,v,x,y) = A{x){d^K{e~'^^,x,y)) B{y) (a^A:(e“^-^, y, x)), 

m = \P\ + |(5|. 

Under the change of coordinates {u, v) {t, r) the function E transforms to H which has the form, 
H{t,T,X,y) = t-^-’^(r(l -r))-^/2-mg-d^(^,y)/(4tr(l-r))p(^^^^^^y)^ 

where 

P{t,T,x,y) G (^“([0,1]2 X 

Taking the order Taylor polynomial of P-^{t,T,x,w) in w we have 

PL{t,T,x,w) = ^ P^{t,T,x)w^. 

\a\<L 

The part of corresponding to Pl is El l, 

(6.17) HL{t,T,X,w) = t — —(r(l-T))-”/2-mg-|«;|V(4tr(l-r)) ^ P^{t, T, x)w^^. 

\a\<L 

Taking the Fourier transform in w, we get the symbol a{HL), 

aiHL,t,r,x,0 = e-l«l'*"(i-")(T(l-T))— P,it,T,x) Y. 

\a\<L FI<ImI 

li^+Zil even 

= ^ RL,^{t,T,x)C- 

H<L 

Now we start the proof of Lemma 6.3 by showing that if L > 3m then 


RLAt,r,x) G C“([0,1] X (0,1] X n). 

This is equivalent to showing that 

(6.18) G C“([0,1] x (0,1] x H x R^). 
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The key here is to write 


(6.19) E^{t,T.x,() 


E 

l/i+^'+t+zc I <m 


{d:,d^K^{e-^^^,x,w)) ,w,x) ), 


where 

Sfx,u,.Ax,w) G ^“([O,!]^ X A), 

and the partial derivatives are taken with respect to the {x, w) coordinates. This is accomplished by 
first changing variables from (x, y) to (x, w), and then using the product rule to take the derivatives off 
the factor ,w,x)). Now from (4.8), 


Hence 

(6.20) E^{t,T,X,i) = i-n-m + IH ^-n/2-m+IH qu^(^ g-|^| V(4ir(l-r)) X^{t,T,X,w) ), 

\u\Kra 

where 

G C°°([0,l]2 X A). 

Denote the order Taylor polynomial of X in re by 

X^{t,T,X,w) = Y 
M<L 

Now replacing X^, in (6.20) by X^ we form 

H^{t,T,X,w) := Y ^-n/2-m+H _ ^yn/2 g_|^|V( 4 tr(l-r)) X^{t,T,X,w) ). 

\u\Kra 

Comparing this with Hl, we find that only terms with degree close to L contribute to the difference 

H^{t,T,x,w) — HL{t,T^x,w). 

Indeed, this difference has the form 

f-"-"^(T(l-T))-’^/ 2 -m g-|^„|V(4ir(l-r)) ^ Yg{t,T,x)w^, 

L — m<\6\<L-\-m 


where 

Yg{t,T,x) G C°°([0,l]2 XD). 
Then computing the symbol by following (6.14), 

a{H^,t,T,x,i) - a{HL,t,T,x,0 
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is a linear combination of terms of the form 

L — m<\6\<L-\-m 

and this is a linear combination of terms of the form 

r, x)^, 

where 

L-m<\e\<L + m, (|0| + |i/|)/2 G N. 

When L > 3m, we conclude that —m + (|0| + |z^|)/2 > 0 and so 

r/2+"*el«l'‘"(i-")(u(F^,t,r,x,0 -a{HL,t,T,x,^)) G ^“([0,1]" x x R"). 

However, again following (6.14), we find that 

a{H^,t,T,x,0 

is a linear combination of terms of the form 

^-n/2-m ^-m X,,^{t,T,x) {i^Y G C“([0,l]2 X 52). 

Hence 

r/2+"*el«l'*^(i-^)u(H'^,t,r,x,0 G C°°([0,1] x (0,1] x 11 x R"). 

which proves (6.18). Interchanging the roles of r and 1 — r by taking the derivatives off the factor 
,x,w)) in (6.19), we prove similarly that 

r/2+’"el«l'‘^(i-^)u(i2L, t,r,rc,0 G C°°([0,1] x [0,1) x H x R^). 

This completes the proof of Lemma 6.3. □ 

Proof of Lemma 6.4- It is clear that since R is smooth on [0,1]^ x 51, 

Q{t,T,x) = [ R{t,T,x) dr 

Jo 

is smooth on [0,1] x [0, oo) x 51. The behavior of Q at T = oo is an application of the method of steepest 
descent. Indeed, 


Q{t,T,x) = f e R{t,T,x) dr + f e R{t,T,x) dr. 

J\t-1/2\<1/A J\t-1/2\>1/A 

The hrst term on the right is smooth for (t, 1/T, x) G [0,1] x [0, oo) x H and has zero Taylor expansion 
at 1/T = 0. For the second term, make the change of variables n = Tr(l — r) to get 


Q{t,T,x) 
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where / G C'°“([0,1] x [0, 3/16] x Q). It is straight forward to show that setting e = 1/T, this is smooth 
for (t, e, x) G [0,1] x [0, oo) x and vanishes at £ = 0. This completes the proof of Lemma 6.4. □ 

Proof of Lemma 6.5. We will need some additional notation. We make Taylor expansions of the terms 
qk{T, x) as T ^ oo: 

OO 

qkiT,x) ~ 

e=o 

Then naturally the terms in these expansions also occur in the expansions of q^{t,x) as t ^ 0, 

OO 

k=0 


We also have a Taylor expansion of qk{T,x) as T ^ 0, 


qkiT,x) ~ ^ 


m =0 


Set 


K-l 


QK{t,T,x) = iQ{t,T,x) - 

V fc=0 / 

so Qk{t, T, x) is smooth up to t = 0, and as T ^ oo we have the Taylor expansion 

OO 

QK{t,T,x) ~ 

where 


e=o 


QK{t,x) = ^^q\t,x) - 


We will start by showing how to obtain the asymptotic expansion of a{s,x,f,) which holds in the L°° 
norm. The symbol estimates will be discussed afterwards. Now a{s,x,^) converges and is analytic in 
s for > n/2 — 1 and we obtain an analytic continuation by using the Taylor series of Q{t,T,x) at 
t = T = 0. Indeed, Taylor expanding Q{t,T,x) at t = 0, 


(6.21) cr{s,x,^) = f P '^/‘^+^QK{t,t\f,\^,x)dt + '^ [ x) dt 

Jt =0 ^_Q Jo 


The first integral converges for > n/2 — 1 — K. Changing variables, 
“1 /"ISI 


/ p-^/^^’^qk{\C\H,x)dt = \^\n-2s-2-2k 

Jo 

Fixing k < 2L, we will analyze 

Ms,x ,0 ■= / T^-^^^+^qkiT,x)dT. 

Jo 


r\i\ 

/ T^-^/^^’^qk{T,x)dT. 

Jo 
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This has an analytic continuation which is seen by making a Taylor expansion at T = 0. Indeed, 


r^-/2+x I qk,m{x)T^j dT 

K-l 

+ E 


Qk,m\^) 


m=0 


s — nj^ + k + m + 1 


lel 


2s—n+2fcH-2m+2 


The integral on the right hand side converges and is analytic in s for JRs > n/2 — K — 1, and the 
other terms are meromorphic with simple poles in the set nj^ — 1 — N"*". This completes the proof that 
has a analytic continuation. Furthermore, set 

poo 

J\i? 

This is convergent for JRs < n/2 — k, and it has an analytic continuation, which can be seen by taking 
a Taylor series as T ^ oo. Indeed, 


(6.22) / r"-^/2+'=gfc(T,x)dT 


L-l 


'i?r 


^s-n/2 + k j ^ j 

qi{x) 


^=0 

L-l 


+ E 


^ s — n /2 + k — i 


ICI 


2s—n-\-2k — 2£ 


i=0 


The first integral is convergent for 3f?s < n/2 — k + L, and is n+ 2 k 2 L^ other terms 

are meromorphic with simple poles in the set n/2 + N. We see that as |.^| ^ oo. 



T^-^/^+'^qk{T,x)dT 




E 

i=0 


qi{x) 


s — n /2 + k — £ 


lei 


2s—n+2fc—2-£ 


and this asymptotic formula still holds after analytic continuation. Now from this we see that 


Us,k{x) = (l)k{s,x,^) + 'il;k{s,x,^) 


pOO 

/ T^-^/^+’^qk{T,x)dT, 

Jo 


which converges for n/2 — /c — 1 < JRs < n/2 — k, has an analytic continuation to s G C with simple 
poles in n/2 + N, and as |^| ^ oo. 


(f>k{s,x,i) 


Us,k{x) + 


qi{x) 


s — n/2 + k — £ 


lei 


2s—n+2fc—2-£ 


£=0 


and from this, 

(6.23) f t^-^/^+\km\x)dt 
Jo 


1^1 


n—2s—2 — 2k 


(l>k{s,x,^) 


1^1 


n —2s —2 —2k 


Us,kix) + 


qHx) 


oo 

^ s-n /2 + /c-. 


ICI 


-2i-2 
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To obtain the asymptotic behavior of a{s,x,^) it remains to compute the asymptotic behavior of 


We deal with this by making a Taylor expansion of QK{t,T,x) as T ^ oo: 


(6.24) f "'^‘^~^^QKit,\^\^t,x)dt 
Jo 






The first integral on the right is provided 3f?s > n/2 — K + L. Choosing K = 2L, we find 

that it is > nl2 — L. Putting this together with (6.23) we find that on the set 

iPs > n/2 — L, as ^ ^ oo, 


~ ICI 


2 L- 1 L-1 


n—2s—2 —2k 


Us,k{x) + E E 


k=o e=o 


s — n/2 + k — £ 


Qidt.x) = -^(i‘{*,x) - , 


We note that 


and so 


(^is,x,0 ~ E 1^1 


2L-1L-1 


n—2s—2 — 2k^ 


s,k{x) + X] E 


k=0 £=0 


s — n/2 + k — £ 


+ Y [ (q^ 




2L-1 \ 

Q\t,x) - Y 

k=0 J 


(6.25) 


Eifi 


n—2s—2 — 2k^ 


s,k{x) + E ^^Mx) ICI ^ + 0(|C| ^) 


where Vg^kix) is extended from its domain of convergence by analytic continuation. It remains to show 
that all the above asymptotic formulas actually hold in the sense of symbols. The key formulas are 
(6.22) and (6.24). For (6.24), we set 


i?^(t,r,x) = QKit,T,x)-YQKit,x)T- 
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is a sum of terms of the form 


^ JO 




where I 7 I = 2j — |/3|. This is bounded by 

|^||7|-2L-2j-2 ^3i.s-n/2+K-L-l 

Jo 

which is equal to 

C(n,s)|^|-2^-2-l/5| 

provided ^R.s > n/2 — K + L. This shows that 




is a symbol of order —2L — 2 when iRs > n/2 — L. 
Similarly, for (6.22), set 


Then 


rkiT,x) = qk{T,x) -'^qi{x)T- 




du / r*-"/"+V^(r,x)dr = -2U\i\^^-^+^\^{\i\\x) 


is a symbol of order 2il?s — n + 2A: — 2L — 1, and 


pOO 

ga / < C'|^|25R*-n+2fc- 

J|£P 


It follows that 


(f>k{s,X,^) - Us^k{x) + 


2s—n-\-2k — 2i _ 


T"-”/2+V^(r,x)dr 


is a symbol of order 23fts — n + 2k — 2L, and putting this together, the 0(|.^| error in (6.25) is 

actually a symbol of order —2L — 2. This completes the proof of Lemma 6.5. □ 
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